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In this paper we describe the physical processes that lead to the generation of Giant Electro- 
Magnetic Pulses (GEMP) on powerful laser facilities. Our study is based on experimental mea¬ 
surements of both the charging of a solid target irradiated by an ultra-short, ultra-intense laser 
and the detection of the electromagnetic emission in the GHz domain. An unambiguous correlation 
between the neutralisation current in the target holder and the electromagnetic emission shows that 
the source of the GEMP is the remaining positive charge inside the target after the escape of fast 
electrons accelerated by the ultra-intense laser. A simple model for calculating this charge in the 
thick target case is presented. From this model and knowing the geometry of the target holder, 
it becomes possible to estimate the intensity and the dominant frequencies of the GEMP on any 
facility. 


INTRODUCTION 

The continuous progress made in the construction of 
high intensity laser opens new domains of physics and 
impacts applications from inertial confinement fusion to 
laboratory astrophysics and material processing. How¬ 
ever the interaction of a high intensity laser pulse with 
a solid target is associated with the generation of in¬ 
tense broadband electromagnetic pulses across a wide 
frequency range from radio frequencies [T] to x-rays [5]. 
With the new generation of PetaWatt lasers, one would 
expect Giant Electro-Magnetic Pulses (GEMP) in the 
Giga- to the Tera-Hertz domain, exceeding 1 MV.m -1 
mm at distances of 1 m from the target, which could be 
very destructive for any electronic device. On the other 
hand, the process responsible for this violent emission, if 
properly controlled, can lead to the production of enor¬ 
mous quasi-static magnetic fields exceeding 1 kT 00 in 
a 1 mm 3 volume, which presents exciting new opportuni¬ 
ties for many applications or fundamental research fields 
such as particle guiding, atomic physics or magnetohy¬ 
drodynamics. 

In this paper, we present the physics of GEMP gener¬ 
ation with short laser pulses and validate it unambigu¬ 
ously in dedicated experiments. In the experimental sec¬ 
tion, we will show that the laser pulse generates a lo¬ 
calised charge Q on the target by the escaped electrons. 
A discharge current I comes from the ground through the 
holder to neutralize it. The GEMP is thereafter emit¬ 
ted by the target holder which behaves like an antenna 
where the discharge current oscillates. This explains two 
major aspects of the GEMP emission. The first is that 
the emission spectrum is defined by the geometry of the 
antenna composed of the target and its support. The 
second is that the strength of the GEMP emission is di¬ 
rectly related to the amount of accumulated charge Q 
in the target. Assuming the target-holder impedanceis 
known, the prediction of this charge becomes crucial for 


experimental designs, either for GEMP mitigation or for 
new conceptions of strong field generation. 

For this reason, in the second part of this article, we will 
present a simplified model to calculate the remaining pos¬ 
itive charge Q in the case of a thick target irradiated by 
an intense laser. We will discuss the domain in which 
our model is valid showing that it may apply to a broad 
range of intense laser facilities worldwide. 

To conclude, we will discuss the opportunity of control¬ 
ling the discharge process to either decrease or modulate 
the GEMP emission, or generate strong quasi-static mag¬ 
netic fields. 


TARGET CHARGING EXPERIMENT 
Experimental set-up 

In our experiment, a Ti:Sapphire laser with a central 
wavelength of 807 nm was focused down to a circular 
Gaussian focal spot of 12 pm FWHM (Full Width at Half 
Maximum) at normal incidence on the flat surface of a 
cylindrical (d = 10 mm) thick (3 mm) target of either 
Copper, Aluminium or Tantalum. The laser temporal 
contrast was 10~ 7 and we varied both the energy on tar¬ 
get [20 mJ—»80 mJ] and the pulse duration [30 fs—>5 ps] 
for each type of target. As shown in Fig [TJ the target 
was supported by a brass wire of 1 mm diameter fixed 
to a connector located in the middle of a large metallic 
ground plane of 20 cm x 20 cm. With such a system 
we could directly connect the brass wire to the core of a 
coaxial cable and fix the distance between the target side 
and the ground plate to l = 4.5 cm. The coaxial cable 
was then connected with a 60 dB attenuation to a 6 GHz 
bandwidth oscilloscope to measure the recharge current / 
circulating in the target. The integration of this current 
provides the charge of the target Q exp generated by the 
pulse. 
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Figure 1. Experimental set-up: the target, its holder and the 
diagnostics. 


A magnetic probe "Bdot", set behind the target at 
a distance of 225 mm from the target center (46 mm 
top and 90 mm right with respect to the laser axis), con¬ 
nected to the same 6 GHz oscilloscope, measures the time 
derivative of the emitted magnetic field in the polarisa¬ 
tion direction perpendicular to the brass wire axis. We 
record B exp , the maximum of the emitted magnetic field 
obtained after integration. 


Experimental results 

We start by focusing our attention on signal spectra. 
Figure [2^ is obtained with a FFT of the raw current sig¬ 
nal and Figure is the FFT of the raw Bdot signal 
integrated in the Fourier space. Despite the various laser 
parameters and target materials used, all electric cur¬ 
rent and magnetic field spectra present the same main 
frequency f exp = 1.0 ± 0.1 GHz. This can be explained 
by considering the target-holder system as a dipole an¬ 
tenna where the ground plate plays the role of a mir¬ 
ror. In that case, the emission frequency is given by: 
f a ~ c/4 {l + ird/2) = 1.2 GHz, where the target-holder 
system is assimilated to a stalk of a length which equals 
the holder length plus half the target perimeter. This es¬ 
timate corresponds fairly well to the measured frequency 
considering the antenna volume simplification. For a bet¬ 
ter definition of the emission spectrum, it is necessary to 
perform a numerical simulation such as those described 
in (7], which reproduces the experimental signal. 

To have a better understanding of the electromagnetic 
emission, we also provide a movie of the electric field 
emitted during the hot electron ejection and the charge 
neutralisation, see Ref. [5]. The first spherical front is 
generated by the charge separation during electron ejec¬ 
tion mm- As the frequency range of this emission is in 
the THz range, it is out of the scope of our diagnostics 
and it is not critical for electronic devices. Note that the 



Figure 2. Typical signal spectra for a) the current and b) 
the magnetic field. f exp = 1.0 ±0.1 GHz, /i == 0.6 GHz and 
j *2 = 1.7 GHz for all experimental points. The black signal is 
the integration noise. 
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Figure 3. Scheme of the electro-magnetic emission from the 
holder. 


dipole formed by the positive charge on the target and 
the ejected electrons is oriented along the laser axis, as 
shown in [5|. The second emission front appears later, 
during the charge neutralisation. The dipole is now de¬ 
fined by the target-holder system and its reflection in 
the ground plate mirror. Figure [3] presents a scheme of 
this emission, while the oscillations of the neutralisation 
current can be observed in the movie [5]. 

We focus now on the amplitude of the electromagnetic 
field. The relation between Q exp and B exp is demon¬ 
strated in Fig. EH+,* , x], for various laser parameters 
and target materials. As the spectra analysis and the 
movie suggest, GEMP are produced by a dipole antenna 
emission. The contact point between the ground plate 
and the brass wire is the center of the antenna and it is 
possible to estimate the amplitude of the magnetic field 
detected by the "Bdot" probe measurement using the far- 
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Figure 4. a) The maximum emitted magnetic field versus the 
associated target charge for several laser parameters and dif¬ 
ferent target materials (copper, aliminum, tantalum), exper¬ 
iments and model (details in text). The experimental points 
are averaged over 5 shots. 


field emission of a dipole antenna at a distance D and at 
an angle 9 with respect to the antenna axis El: 


B 


a 


p 0 I cos (f cos 0) 
27 tD sin 9 


(1) 


In our case D = 240 mm, 9 = 77° and the discharge cur¬ 
rent is estimated by I = 1.6 cQ exp / ( l + ttcL/2) where 1.6 
is the impedance adaptation coefficient. As the "Bdot" 
probe signal is noisy (in black in Fig. §>) and also con¬ 
tains cavity modes of the experimental chamber (in blue 
in Fig. §>), we filtered it between fi and / 2 to isolate 
the antenna component signal (in green in Fig. ©>)■ 
The experimental filtered B exp ([o] in Fig. |4j) is com¬ 
pared to the magnetic field from Eq. ([l]) ([solid line] in 
Fig. The closed agreement shows that the ground- 
holder-target system is the antenna, which emits GEMP 
and that the accumulated charge is the energy source of 
GEMP. Thus, a model was developped to estimate this 
accumulated charge. 


TARGET CHARGING MODEL 
Model assumptions 

The evaluation of the charge Q accumulated on the 
target after the end of a laser pulse is a multi-physical 
problem, which involves various effects related to laser- 
matter interaction, particle diffusion and collision in the 
target, electromagnetic properties and an overall descrip¬ 
tion of the chamber environment [zma. Here, we present 
a model, which estimates, to the nearest order of mag¬ 
nitude, the target charge assuming, that the laser ab¬ 


sorption coefficient is known. It includes all major phys¬ 
ical effects and is applied to thick metallic targets, i.e. 
the target thickness is larger than the range of acceler¬ 
ated electrons. Thus the electrons are ejected only from 
surface irradiated by the laser. The basic principles of 
the electron dynamics can be described as follows. The 
laser pulse accelerates electrons at the impacted target 
surface. They spread into the target dissipating their 
energy inside and creating an electric potential at the 
target surface H3HI3- This potential is a barrier that 
the accelerated electrons must overcome in order to es¬ 
cape completely from the target. The number of escaped 
electrons defines the net positive charge left on the target 
surface. 

Instead of following the dynamics of the hot electrons, 
as was done in Ref. [7j, here we focus our attention on the 
maximum number of electrons that may escape. The ba¬ 
sic assumptions concern the characteristics of the accel¬ 
erated electrons. We assume that they are created in the 
laser focal spot of radius ri as , and their density equals ap¬ 
proximately the relativistically corrected critical density 
n c = l.l^/l + ao 2 /2 ■ 10 21 Aj“g cm -3 , with Ai as being the 
laser wavelength in //in, a o = O.SSAias/j,),, the dimension¬ 
less laser vector potential and /i as = 0 . 65 -Ei as /ti as 7 rrj 2 as 
the laser intensity with units of 10 18 WcnG 2 [T6j. The 
energy of the accelerated electrons follows a Maxwell- 
Jiittner distribution function ©> characterised by the hot 
electron temperature T^: 

fe{e e ) = -^ 7 e p e exp(-e e /T /l ), (2) 

where p e is the electron momentum, y e the relativistic 
factor and A the normalisation factor. Such an approx¬ 
imation for the electron energy distribution is reason¬ 
able because the laser accelerated electrons are scattered 
by the electrons and ions in the target. For the tem¬ 
perature, we use Beg’s empirical law [ T8] in the interval 
0.03 < a o < 1 or ponderomotive scaling m for higher 
laser intensities: 

Th = m e c 2 max |o. 47 a^ 3 , \J 1 + — 1 j . (3) 

For intensities below 10 15 Wcm~ 2 (ao < 0.03 ), the hot 
electron temperature is estimated from the model of laser 
collisional absorption [2(J|, T h = 3m e c 2 aQ^ 3 . The ratio 
between the absorbed laser pulse energy ?/ a b s -Ejas and the 
mean energy of accelerated electrons (e e ) = f Q £ e f e de e 
gives the total number of accelerated electrons: 

-^e = Aibs Alas / (a) (4) 


Potential barrier calculation 

The critical section of the model is the estimation of 
the potential barrier, which eventually defines the max- 
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inrum number of escaped electrons. The first contribu¬ 
tion to this barrier comes from the cloud of energetic 
electrons outside that target, which creates a charge sep¬ 
aration potential (f> t h as described in the ion acceleration 
model HMg. The high energy electrons may overcome 
the potential barrier and escape completely from the tar¬ 
get. The remaining positive charge generates a global 
background potential 4>bg > which spreads over the target 
surface. First, we focus on the first part of the poten¬ 
tial, (f> t h, which is described by the Poisson equation © 
assuming a Boltzmann distribution of electrons in the 
potential barrier: 

e 0 A (t) th = - e(rii - n c exp(e^ tft ,/T/ 1 )), (5) 


where e is the electron charge and eo is the vacuum di¬ 
electric permittivity. The ion density is described by a 
Heaviside function, n, = n c H(x), where x = 0 defines 
the target surface. This hypothesis is valid as long as the 
characteristic scale of the ion density is smaller than the 
hot electron Debye length defined below. 

This equation has a divergent solution for the poten¬ 
tial in one dimension eh. Here, we present a method 
which determines completely the potential barrier: (j) t h 
is convergent inside and outside the target. In contrast 
to the ion acceleration problem EMU, where only the 
electric field at the surface matters, the problem of elec¬ 
tron escape requires a knowledge of the whole potential 
profile. In the one-dimensional model, one obtains the 
following electron density distribution: 


n c exp ( — exp(— ^x/Xd)) x > 0, 
n c (e 1 / 2 — x/\dV%) 2 x < 0 


where Xjj = \JeoTh/n c e 2 is the hot electron Debye 
length. The coefficient £ = 0.9288 is calculated to con¬ 
serve the electro-neutrality. As the electron density © 
is a converging function: we can use it to determine the 
potential (j>th with a three-dimensional calculation. The 
hot electrons are distributed over a cylinder of radius 7'i as 
with the i-axis perpendicular to the target surface (see 
figure [5j . The potential is numerically calculated from 
the charge distribution Eq. © : 



Figure 5. Example of a space charge reconstruction for the 
normalized radius R = 50. Top inset: axial distribution of 
the ion and electron density from Eq. ((|6|). Bottom inset: 
the axial potential. 



4>th(r,x) 


e 

47re 0 


( rii(x ') — n e (x')) dx' r' dr' dd 



r'r cos 6 


. (?) 

It is important to note here that the potential normalized 
to the hot electron temperature, (f> t h = e^th/Th, depends 
only on the normalized radius R = r\ as /\D- Once the po¬ 
tential (f>th (f, 0, x) is calculated, we extract the normal¬ 
ized potential barrier labelled AIt is defined as the 
maximum potential variation over the x-axis, averaged 
over the cylinder section. For large values of R > 10, we 
provide a fit for the potential barrier at the center of the 


Figure 6. Dependence of the maximum normalized potential 
barrier A <j>th on the normalized laser spot radius R. 


target cylinder 

A %th ^ 41og fl (1 - \og R /(I + logfl)) (8) 

with log^j = log 10 R. The potential barriers are presented 
in Fig. [S] 
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Charge estimations 

We now combine the various aspects of the model to 
estimate the target charge after the laser pulse. Knowing 
the temperature ©> , the normalized potential barrier <©> 
the distribution function ([2]), and the total number of 
hot electron Q, one can readily determine the amount 
of charge escaping the potential barrier: 

/»oo 

Q = eN e ^ f e (£ e )d£ e . (9) 

j ThAffith. 

The model’s assumptions are valid for several laser fa¬ 
cilities with various laser parameters, with the exception 
of high contrast systems where the electron acceleration 
may not be descripted by ^. 


Effect of the target size on the charge 


We focus now on the background potential (f>b g , which 
is due to the uncompensated positive charge left at the 
target surface by the escaped electrons. This positive 
charge spreads over the whole target surface, which re¬ 
veals a target size effect. For infinite or big enough tar¬ 
gets, the background potential remains negligible because 
the surfacic charge is low. For small targets, the charge 
concentrates at the target surface and develops a back¬ 
ground potential, which may overcome the potential bar¬ 
rier. To investigate the impact of the background poten¬ 
tial on the target charge, we estimate 4>b g as follows: 


4>b g — 


Q 

27re 0 d 


( 10 ) 


This formula assumes a uniformly charged disc of a neg¬ 
ligible thickness compared to its radius. In the domain 
of practical interest, the disc radius is much larger than 
the hot electron localisation zone, which implies that <f>b g 
is fairly constant along the laser axis in the <f>th barrier 
range. As a consequence, the whole potential curve tfith 
is offset by the background potential. Two examples of 
the total normalized potential (f> t h + 4’b g are shown in 
Fig. 0 If (f>b g < (f>mini the potential barrier value is 
unchanged (red dashed line). If (f>b g > (limin', the poten¬ 
tial barrier increases (blue solid line). As shown in ([9]), 
the potential value is the limit of integration of the hot 
electron distribution function and therefor the chatge is 
highly dependent on this parameter. In fact, the increase 
in the potential barrier inhibits the ejection of electrons 
and limits the charging process. This strangled regime 
appears for small targets as they concentrate the sur¬ 
facic charge and enhance the background potential. This 
strangled regime also applies for low laser intensities be¬ 
cause of a low value of <t> m in ■ As the ejection current 
strangling appears when (f>b g = 4> m in , we use this rela¬ 
tion to define the target size where the strangling effect 



Figure 7. The potential for R = 50 with an effective (solid 
blue line) or ineffective (dashed red line) arbitrary background 
potential d = d/\D = 1000. 


is suppressed: 

, _ Q 

Q"min — 0 , 

Z7T6o(pmin 


(ii) 


If the target is too small, we can also make a charging 
prediction, assuming an instantaneous struggling: 


Qbg 2lTdeo(f> min- 


( 12 ) 


In our experiments, the largest critical size is d m i n = 
5 mm, which is still smaller than the target size. While 
the strangled regime does not apply to our experiment, it 
was observed in Ref. [ 23 ] and in Ref. [ 21 ] where d m i n was 
larger than the target diameter of 1 cm. Consequently, 
the authors measured a GEMP amplitude proportional 
to the target diameter, in agreement with Eq. (12). 


COMPARISON BETWEEN MODEL AND 
EXPERIMENTS 

A comparison between the measurements and the 
model predictions is shown in Fig. s* for scans of laser 
energy E \ as and in Fig. [8 ]d for scans of pulse duration 
fi as . The first important result is that the charge is inde¬ 
pendent of the target material. This feature is especially 
obvious at low laser energy and is implicit in the model, 
which does not include any target parameters (with the 
exception of the absorption coefficient). The second re¬ 
sult is the relative independence of the charge with re¬ 
spect to the laser pulse duration. A physical explana¬ 
tion of this independence can be seen undertood if we 
focus on the number of accelerated electrons given by 
©: for short pulses, the escaping electrons represent a 
large fraction of the small total number N e of hot elec¬ 
trons, while for longer pulses, they are a small fraction 
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Figure 8. Target charge in nC measured experimentally (dots) 
and calculated with the model (line) as a function of laser 
energy £) as (a) and pulse duration ti as (b), for different target 
materials. The points are averaged over 5 shots. The error 
bars are the standard deviation on these points. 


of the numerous hot electrons. Also, it should be noted 
that the charge scales linearly with i?i as as in f22j. The 
model’s predictions are in reasonable agreement with the 
experimental data for a broad range of laser parameters, 
considering the simplicity of the model. Nevertheless, 
there is some behaviour linked to the target material, for 
pulses longer than 0.5 ps, as illustrated in Fig. [8]d. As 
the material affects the absorption coefficient, this de¬ 
pendence may be accounted for by an appropriate vari¬ 
ation of the laser energy absorption, taking into account 
a broad range of pulse durations. 

While the model is only corroborated experimentally 
for low laser energies, the model’s assumptions are valid 
in a broader laser parameter space. Thus we extend the 
predictions of accumulated charge Q as a function of laser 
energy and pulse duration in Fig. [9] for a laser focal 
spot radius of 6 pm and an absorption coefficient of 40%. 
These results should be carefully considered at high laser 
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Figure 9. Expected target charge Qt in nC and in the loga¬ 
rithmic scale calculated from the model in function of the laser 
energy and the pulse duration for ri aa = 6 pm, r; a b 3 = 0.4, and 
Aiaa = 0.8 pm. The dashed line shows the domain explored in 
the present experiment. 


energy far from the validated range. This map highlights 
two different regimes of target charging. First, an almost 
complete hot electron ejection takes place if (e e ) > A (f> t h- 
In this case the target charge can be approximated by 
Q ~ eN e . Second, there is a quasi-stationary regime 
where the laser pulse duration is longer than the hot elec¬ 
tron cooling time |7]. In this case, the current of ejected 
electrons is equal to I e = Q/ti as . 


CONCLUSION 


In conclusion, by comparing our analytical model with 
our experimental results, we demonstrate that GEMP 
emission proceeds in two steps: the target is charged by 
the laser and is subsequently discharged through the tar¬ 
get holder. 

In consequence, an ideally isolated target cannot generate 
GEMP in the GHz domain, even though it accumulates 
a charge. It may, however, produce TeraHertz waves due 
to the ejected electrons, x-rays due to the target plasma, 
or weak GEMP due to the polarisation current in the di¬ 
electric holder. Our second conclusion is that this simple 
fact allows us to determine two kind of laser-matter inter¬ 
action experiments. If the charging time is shorter than 
the neutralisation time, there is a charge accumulation 
and then a GEMP emission. When these characteristic 
time are inversed, there is no charge accumulation but 
rather a constant current, which generates weak GEMP. 
Assuming a neutralisation time of l/c, a charging time 
proportional to ti as and the holder size around 10 cm, 
laser pulses shorter than 300 ps produce GEMP with 
a strength proportional to the laser energy, while laser 
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pulses longer than 300 ps generate GEMP strongly in¬ 
hibited by the neutralisation current. Our experiment 
produced GEMP efficiently because the ps and sub-ps 
laser pulses were much shorter than the characteristic 
discharge time (around 0.1 ns). Moreover, by properly 
shaping the target holder, it is possible to control the 
GEMP spectrum for diagnostics protection or for gener¬ 
ation of strong quasi-static magnetic fields M. 
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